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ABSTRACT 

We present the results of a series of calculations studying the collapse of molecular cloud 
cores performed using a three-dimensional smoothed particle hydrodynamics code with ra- 
diative transfer in the flux-limited diffusion approximation. The opacities and specific heat 
capacities are identical for each calculation. However, we find that the temperature evolu- 
tion during the simulations varies significantly when starting from different initial conditions. 
Even spherically-symmetric clouds with different initial densities show markedly different 
development. We conclude that simple barotropic equations of state like those used in some 
previous calculations provide at best a crude approximation to the thermal behaviour of the 
gas. Radiative transfer is necessary to obtain accurate temperatures. 
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1 INTRODUCTION 

Radiative transfer is an important phenomena in star formation. 
Radiation sets the temperature of the gas during the collapse of a 
molecular cloud core. This both influences the degree of fragmenta- 
tion of the cloud, and sets the minimum mass of brown dwarfs (the 
opacity limit for fragmentation: Ilow & Lvnden-Belllll97rjt) . Once 
protostars have formed in a cloud, radiative and mechanical feed- 
back from them can affect subsequent star formation. Such feed- 
back mechanisms include protostellar jets and outflows from low- 
mass stars, and ionisation from massive stars which creates HII re- 
gions and destroys a cloud. 

Computer simulations are vital in our efforts to understand 
the complex problem of star formation. Many previous simula- 
tions have used the smoothed particle hydrodynamics method 
(SPH) (e.g. " 



Bonnell et al 



sed the smoothed particle hydrodynamics n 
Gh^ld^Mona^harj |l98l|: [PoMracketalJ 
Il99lt iNelson & Papaloizo3ll993t iBate et all 
F" 



Klessen et al 



1991 



1995 



1998). Other methods used are typically based 
aroun d grid-based code s (e.g. [Larson 1969; Boss & Bodenheimer 
[ 19791: iBoss & Mvhilll [l99a iBurkert & Bodenheimed Il993t 
bruelove et alJ Il997t Boss et alj l200(j) . SPH is a Lagrangian 
method first developed bvlLTc^TIl977l) and fGingold & Monaehanl 
h977l) (see iMonaghanl fl992L for a review). It approximates the 
fluid as a series of discrete fluid elements denoted by individual 
SPH particles and uses interpolation to obtain the fluid variables 
at any point in the simulation. SPH is conceptually simple to 
understand, and can naturally adapt its resolution to the local 
density distribution, unlike grid-based codes which require com- 
plex adaptive-mesh refinement algorithms to perform the same 
task. This property makes it ideal for use in star formation, where 
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densities may range over many orders of magnitude in a single 
simulation. 

Despite these advantages, few attempts have been mad e to in- 
clude radiative transfer into SPH l Luc vll 19771 lBrookshawlll985l 
Il98dl0xlev & Woolfsonl2003tlWhitehouse & Batel2004 . and un- 
til recently iBastien. Cha & Viadl2004 120051) SPH with radiative 
transfer has not been applied to star formation. Instead, many past 
simulations have simply used isothermal or barotropic equations of 
state to model the collapse of a molecular cloud. The former is only 
valid up to densities of ~ 1CT 13 g cm" 3 at which point the cloud 
traps radiation efficiently enough for the cloud to begin to heat 
up. The latter is usually based on the evolution of the temperature 
at the highest density during the collapse of spheri cally symmet- 
ric clouds as calculated using radiative transfer (e.g. | Larsonll969t 
Winkler & Newman 1980; Masunaga & Inutsuka 2000). However, 
a barotropic equation of state can at best only hope to provide an 
adequate description of temperature at the density maximum; it is 
unlikely to give an accurate temperature distribution during a three- 
dimensional calculation with c omplex density and velocity struc- 
ture. Indeed, Boss et al. (2000) performed grid-based calculations 
of the collapse of a molecular cloud core both with a barotropic 
equation of state and with radiative transfer in the Eddington ap- 
proximation and found they differed somewhat. However, they did 
not examine in detail how the relation between temperature and 
density differed from that of the barotropic equation of state spa- 
tially and temporally, or its dependence on initial conditions. 

Whitehouse, Bate & Monaghan 120051) recently presented an 
implicit algorithm for calculating radiative transfer using the flux- 
limited diffusion approximation within the SPH formalism. This 
paper describes a three-dimensional implementation of this algo- 
rithm and uses it to examine the thermodynamics during the col- 
lapse of molecular cloud cores. Section [2| describes the changes 
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necess ary to the radiative transfer algorithm of [Whitehouse et all 
(2005) for use in three dimensions and the initial conditions for 
our star formation calculations. Section [3] presents the results of 
simulations of the collapse of molecular cloud cores with different 
initial conditions and examines the evolution of their temperature 
structure. Finally, section summarises the main conclusions of 
this paper. 



2 METHOD AND INITIAL CONDITIONS 

The code used in this paper is based on that of Bate 



Bate et al. 1995), which originated from that of Benz l 
Benz et alll 990). The code includes individual timesteps, and uses 
a tree to calculate self-gravity. The smoothing lengths of particles 
are variable in time and space, subject to the constraint that the 
number of neighbours for each particle must remain approximately 
constant at N ne \ g <n = 50. We use the standard form of artificial vis- 
cosity (Monaghan & Gingold 1983; Monaghan 1992) with strength 
parameters a w = 1 and B v = 2. To perfor m radiation transport, the 
algorithm from Whitehous e et alJ 12005T) . adapted for use in three 
dimensions, was added. The code has been parallelised by M. Bate 
using OpenMP. 



2.1 Three-dimensional flux-limited diffusion 

In a frame co-moving with the fluid, and assuming local thermal 
equilibrium (LTE), the equations governing the time-evolution of 
radiation hydrodynamics (RHD), integrated over frequency, are 
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(Mihalas & Mihalasl fl98l iTurner & Stonell200ir) . In these equa- 
tions, D/Dt = d/dt + v • V is the convective derivative. The sym- 
bols p, e, v and p represent the material mass density, energy den- 
sity, velocity, and scalar isotropic pressure respectively. The total 
frequency-integrated radiation energy density, momentum density 
(flux) and pressure tensor are represented by E, F, and P, respec- 
tively. In this paper, we use a grey opacity, k (i.e., it is independent 
of frequency). 

To solve radiation transport within SPH, we evolve both the 
specific internal energy of the gas u, and the specific radiation en- 
ergy £ = Ejp implicitly using the algorithm of Whi tehouse et alJ 
(2005). The principle difference between the one-dimensional and 
three-dimensional forms of the radiative transfer equations is the 
radiation pressure term (the second term on the right-hand side of 
equation [3}. In one-dimension this term can be written using the 
divergence of the gas velocity. However, in three dimensions this 
becomes the tensor product Vv:P, where the Eddington pressure 
tensor P is given by 



The components of the Eddington tensor / are given by 

/= 1(1 -/)/+! (3/- (7) 

Here / is the dimensionless scalar Eddington factor, and h is the 
unit vector in the direction of the energy density gradient VE/\VE\, 
For example, in Cartesian coordinates, the first component of the 
tensor is given by 

1 2(3/-D(f) 2 
2 ( f> + 



|V£p 



(8) 



and the next by 



(3/-D(f)(f) 



(9) 



|V£p 
and so on. 

The scalar Eddington factor / is related to the flux-limiter A, 
by the expression 

f = A + A 2 R 2 , (10) 

where/? = \VE\/(kpE). In this p aper, we choose the flux limiter of 
lLevermore"& Pomraninal 1 198 ll> 

2 + R 
6 + 3R + R 2 

In the optically thick limit, R — > 0, A — > 1/3, and P becomes 
isotropic. 

2.2 Specific heat capacity and opacity 

The above equations are closed by the application of an equation of 
state for the gas. We use the ideal gas equation of state 

P = —pT e , (12) 
P 

where R g is the gas constant, p is the mean molecular mass and the 
gas temperature is T g = u/c\. The temperature of the radiation is 
given by T T = (E/a) l/4 . 

W e use the specific heat capacity c v of lBlack & Bodenheimed 
( 1975), which accounts for the dissociation of molecular hydrogen, 
and the ionization of both hydrogen and helium. It omits any con- 
tribution due to metals. 

c v = X(l-y)E(H 2 ) + [1.5X(l+x)y + 0.375Y(l+zi+z 1 Z2)]R B 
+X(l.304 x 10 13 x + 2.143 x 10 12 )y/r g + 

y (5.888 x 10 12 (1 -z 2 )+ 1-892 x 10 13 z 2 ) Zl /r g , (13) 

where X and Y are the mass fractions of hydrogen and helium re- 
spectively (in the simulations presented in this paper, X = 0.70 
and Y = 0.28), y is the dissociation fraction of hydrogen, x the 
ionisation fraction of hydrogen, and z\ and z 2 are the degrees of 
single and double ionisation of helium, respectively. E(H2) gives 
the con tribution to the specific heat capacity from molecular hy- 
drogen <Black & BodenheimeJll975l) . The ionisation fractions are 
calculated using the Saha equation. 

The variati on in mean m olecular mass n with temperature was 
also taken from Blac k & Bodenheimed J 19751) to be 

p- [ = [2X(l + y + 2xy) + Y (1 + Zi + ZiZ 2 )] /4. 



P = fE. 



(6) 



(14) 

The opacity table from Alexander] 1 19751) (the 
fourth King model) was use d for the gas opacity and 
IPollack. McKay & Christoffersorl jl985t) for the dust opacity. 
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The opacity values were stored in a table containing the opacity at 
various temperatures and densities up to 10,000 K and 1 g cm 4 
respectively, and bilinear interpolation in log-space was used to 
get the required value. Above 10,000 K, the opacity was taken to 
be the lesser of extrapolating from the last two points in the table, 
or Kramer's opacity 



k k = 1-2512 x 10 /37"!, 



(15) 



until the electron scattering opacity k cs = 0.4 cm 2 g" 1 becomes the 
dominating opacity source at very high temperatures. 

For temperatures above a few million K, the opacity, specific 
heat capacity and mean molecular mass become constant because 
the contribution of metals is neglected. These three quantities were 
updated during every iteration of the implicit scheme. 



2.3 Initial conditions 

The simulations described in this paper used three different sets of 
initial conditions. The first was the spherically symmetric collapse 
o flBoss & Mvhilll <1992h . They began with a 1 .2 M e sphere of uni- 
form density p = 1.7 x 10"' 9 g cm 4 , giving a free-fall time of 
tff = 1.61 x 10 5 years. The cloud radius was 1.5 x 10 17 cm, and the 
initial temperature of both the gas and the radiation field was 10 K. 

The other two initial conditions were based on those of 
iBoss & BodenheimeJ ll979l) . They began with a 1.0 M Q sphere of 
density p = 1.44 x 10" 17 g cm 4 and radius 3.2 x 10 16 cm, giving 
a free-fall time of t s = 1.75 x 10 4 years. Their cloud was initially 
in solid-body rotation with an angular velocity of 1.6 x 10" 12 rad 
s" 1 . Superimposed on the underlying density was an m = 2 density 
perturbation satisfying 



p = po[l +0.5 cos(mt/>)] 



(16) 



where (f> is the angle in the plane perpendicular to the axis of rota- 
tion. The ratio of thermal energy to magnitude of the gravitational 
energy was initially 0.26 and the ratio of rotational energy to the 
magnitude of the gravitational energy was 0.20. Using our equa- 
tion of state, the initial temperature of the gas and radiation were 
therefore set to 12 K. We performed calculat ions with i nitial condi- 
tions identical to those of lBoss & BodenheimeJ ll979l) . and also a 
spherically symmetric version in which there was no in = 2 density 
perturbation applied and the cloud was not rotating. 

We performed simulations with 5000, 50,000, 150,000 and 
500,000 particles. Of these, simulations with 50,000 particles and 
gr eater resolve the Jeans mass according to the criteria set out 
in iBate & Burkertl ll997l) . The calculations were performed on 
the United Kingdom Astrophysical Fluids Facility (UKAFF). The 
highest resolution (500,000 particle) calculation took a total of ap- 
proximately 3000 CPU hours (running across multiple processors). 



3 RESULTS 

A brief summary of the expected beha viour during a spherically 
symmetric collapse is as follows (see Boden heimer & Sweieartl 
Il968t lLarsortll969l for more details). A Jeans-unstable uniform 
cloud of molecular gas will begin to collapse under its own self- 
gravity. The collapse occurs such that the density in the outer parts 
of the cloud falls off with radius r as p cc r" 2 , while the uniform- 
density inner part of the cloud collapses to higher and higher den- 
sities but contains a decreasing fraction of the mass as the collapse 
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Figure 1. We show the evolution of the maximum temperature and 
density during the ca lculations that start from the initial conditions of 
iBoss&Mvhill 119921) for different resolutions. The calculations were per- 
formed using 5000 (dot-dashed green line), 50,000 (long-dashed blue line), 
150,000 (dotted red line) and 500,000 (solid magenta line) particles. The 
results converge as the resolution is increased, with 50,000 or greater num- 
bers of particles required to achieve reasonable results. For comparison, the 
short-dashed line gives the tem perature-density relation of the barotropic 
equation of state used by Bate Ul993) . 



proceeds. Initially, the collapse is isothermal, as the material is op- 
tically thin. However, once the central region reaches a critical den- 
sity (~ 10~ 13 g cm" 3 ) it becomes optically thick and starts to heat 
as energy can no longer be radiated away. The central temperature 
and density both rise rapidly, until thermal pressure can counter the 
collapse. This results in the formation of a pressure-supported core 
in the centre of the cloud. Mass continues to infall onto this core, 
increasing its mass and temperature. While at low temperatures, 
the molecular hydrogen behaves like a monatomic gas (the ratio of 
specific heats y = 5/3). However, when the temperature approaches 
~ 100 K (at densities ~ 10" 11 g cm" 3 ) additional degrees of free- 
dom become available and y = 7/5. This phenomenon which is in- 
cluded in our specific heat capacity is not typically included in the 
simple barotropic equation of state. Once the core reaches ~ 2000 
K, the molecular hydrogen begins to dissociate. As energy is di- 
verted into this dissociation rather than thermal support, the core 
begins a second collapse phase. Once the hydrogen has been fully 
dissociated, it forms a second ("stellar") core, once again supported 
by thermal pressure. This core continues to increase mass and tem- 
perature as material falls onto it from the envelope. 

We present results showing the evolution with time of our cal- 
culations in FiguresQto[8] The short-dashed bl ack line in th ese fig- 
ures is the barotropic equation of state used bv lBatd 11998). given 
for comparison. Unless otherwise stated, we plot the temperature 
of the gas rather than that of the radiation field. 

As far as we are aware, these are the first three-dimensional 
radiative transfer calculations to follow the collapse of a molecu- 
lar cloud core beyond the formation of the first pressure-supported 
core and the dissociation of molecular hydrogen to the formation of 
the stellar core. lBos'sl fl984) performed one- and two-dimensional 
radiative transfer calculations that followed the collapse of a molec- 
ular cloud core well beyond the formation of the first pressure- 
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Figure 3. Radial profiles of density (left panel) and gas temperature (right panel) at various stages during the highest resolution (500,000 particles) 
iBoss & MvhiJ 1 1992ft calculation. A profile is plotted every order of magnitude in maximum density from 10~ 12 to 10~ 3 g cnT 3 and finally at 2 X 10~ 3 
g cm . For comparison, the thick solid black lines give the profiles from lBoss & Mvhilj <1992l) when the temperature in their calculation was ~ 1000 K. 
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Figure 2. We show the evolution of the maximum gas (solid green line) 
and radiation (dotted red line) temperatures at the maximum density dur- 
ing the lBoss & MvhiJ (l992) calculation. This was our highest resolution 
(500,000 particles) calculation. The gas and radiation temperatures are iden- 
tical except during the transition from optically thin to optically thick. For 
comparison, the short-dashed line gives the temperature-density relation of 
the barotropic equation of state used bv B ateUl998l) . 



supported core. However, in order to accomplish this he had to 
employ several numerical artifices, begging the question of how 
we are able to perform three-dimensional calculations. Boss' first 
artifice was to implement two regions of his grid that were evolved 



using different timesteps in order to follow the first core over many 
dynamical times using short timesteps but simultaeously model the 
envelope with larger timesteps. He also damped out oscillations of 
the first core to decrease the amount of computational effort re- 
quired. In our calculations, we encountered no difficulties in mod- 
elling the cloud collapse through the first core phase and onto the 
formation of the stellar core. We believe there are th ree reasons 
for this. First, o ur code employs individual timesteps iBatdl 19951 
Bate et al. 1995) for each particle (i.e. similar to, but even more ef- 
ficient than, the way Boss evolved his calculation in with two spa- 
tially distinct timesteps). By the end of our highest resolution cal- 
culation, some particles within the stellar core were being evolved 
using timesteps of less than 1/40 of a second (« 10 5 times smaller 
than those in the outer parts of the cloud). Second, our calculations 
used the standard form of SPH viscosity with none of the possible 
viscosity-reducing formulations. In each calculation, the first core 
undergoes oscillations after its formation, as observed in earlier 
work, but the viscosity likely damped these oscillations in a similar 
way to Boss' artifice. Finally, in the twenty years between Boss' 
and our calculations, computers have become very much quicker. 



3.1 Results using Boss and My hill initial conditions 

Figure Q shows the evolution of maximum temperature and den- 
sity during the Boss & Mvhilll il 992h collapse calculations. We per- 
formed the collapse with four different resolutions. With 5000 par- 
ticles (green line) the cloud heats at an earlier stage than in the 
higher resolution simulations, presumably due to insufficient res- 
olution. The 50,000 particle collapse (blue) is much cooler for a 
given maximum density, while the two highest resolution simula- 
tions appear to be converging towards a single curve. We conclude 
that 50,000 to 150,000 particles (i.e. approximately the same num- 
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Figure 4. Snapshots of temperature versus densi ty at various stag es during 
the highest resolution (500,000 particles) Boss & Myhilj 1 1993) collapse. 
A snapshot is plotted every order of magnitude in maximum density from 
10~ 12 to 10~ 8 g cnT 3 and finally at 2 X 10~ 3 g cm 4 . Note that at late times 
the temperature differs by more than an order of magnitude from that given 
by the barotropic equation of state (dashed line). 



ber as required to resolve the Jeans mass) are sufficient to model the 
thermal behaviour reasonably accurately. The evolution of maxi- 
mum temperature and density follows the barotropic equation of 
state (dashed line) in a qualitative sense, but there are still differ- 
ences in temperature of up to « 50 percent between the radiative 
transfer and the barotropic equation of state at various stages dur- 
ing the collapse. 

Figure [2] shows the evolution of both the gas temperature 
(green) and the radiation temperature (red) for the highest reso- 
lution (500,000 particles) IBoss & Mvhiill 1 19921) collapse. The ra- 
diation temperature is equal to the gas temperature at most stages, 
except in the vicinity of the transition between the optically thin 
and thick regimes. This can be seen in the figure from densities of 
~ 10~ 15 to 10~ 12 g cm" 5 when the cloud initially begins to heat. 
Here the radiation temperature is up to a factor of two greater than 
that of the gas. Later in the calculation, gas that is optically thick 
has the same temperature as the radiation and the two remained 
coupled as the gas collapses to very high densities. However, there 
is always a region where the gas transitions from optically thin to 
optically thick where the temperatures differ. 

Figure0shows density (left) and gas temperature (right) plot- 
ted against radius as the collapse passes various values of the max- 
imum density. The first core can be seen at a radius of ~ 10 14 cm 
where the density and temperature gradients change abruptly. Sim- 
ilarly, the transit ion to the stellar core can be seen at a radius of 
~ 10 12 cm. lBoss & Mvhill 1 1992) give the density and temperature 
profiles when the central temperature reaches » 1000 K. We repro- 
duce their profiles as the thick solid black lines in Figure [5] Our 
results are in good agreement with theirs, the main difference be- 
ing that our results give slightly higher temperatures at radii ~ 10 15 
cm, presumably due to the fact that we use flux-limited diffusion 
while lBoss & Mvhilll 1 19921) use the Eddington approximation (the 
former retards radiation transport near the surface of the first core). 

In Figure we plot temperature versus density at a series 
of snapshots during the highest resolution calculation (since both 



Figure 5. The evolution the maximum temperature and density during cal- 
culations starting from three different sets of initial conditions. The three 
calculations are the sphericall y-symmetric (dotted green lin e) and m = 2 
rotating (long-dashed red line) Boss & Bodenheimer 1 1979) collapses and 
the spherically-symmetric Boss & Mvhill 1 1992) collapse (solid blue line). 
All calculations were performed using 50,000 particles. Note that different 
initial conditions can lead to very different temperature evolution. For com- 
parison, the short-dashed line gives the tempera ture-density relation of the 
barotropic equation of state used bv lBateHl998l) . 



temperature and density are monotonically decreasing functions of 
radius, we can plot these as lines for each snapshot). The figure 
clearly shows that even if it were possible to construct a barotropic 
equation of state to follow the evolution of maximum density and 
temperature, it would generally underestimate the temperature of 
the gas away from the maximum. In the figure, the barotropic equa- 
tion underestimates the temperature by more than an order of mag- 
nitude at densities ~ 10~ 13 g crrT 3 at late times. 



3.2 Results using spherically symmetric Boss and 
Bodenheimer initial conditions 

Figure [5] shows the evolution of the maximum value of tem- 
perature and density for the three types of initial conditions, 
each performed with 50,000 particles. The spherically-symmetric 
IBoss & BodenheimeJ 1 19791) coll apse (green line ) begin s to heat 
at a much lower density than the Boss & Mvhill 1 1992) collapse 
(blue). Essentially the only difference between these calculations is 
the initial density of the two clouds, the former is initially almost 
two orders of magnitude denser, and yet the temper ature evol utions 
are qu ite different. Figure|S|shows that, as with the lBoss &l vlvhill 
1 1992) collapse, the spatial and temporal evolution of the tem- 
perature is also complex. Again, these results demonstrate that a 
barotropic equation of state cannot accurately describe the temper- 
ature distribution during such cloud collapses. 



3.3 Binary star formation 

The result of the lBoss & BodenheimeJ J 19791) initial conditions, in- 
cluding the m = 2 density perturbation and rotation, is to form a 
wide binary with a separation of « 1000 AU iBate et alll995l) . Us- 
ing these initial conditions, the evolution of the maximum temper- 
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Figure 6. Snapshots of temperature versus density at various stages dur- 
ing the 50,000 particle spherically-symmetric iBoss &B odenheimerl Il979l) 
calculation. A snapshot is plotted every order of magnitude in maximum 
densit y from 10~ 13 to 10~ 9 g cm -3 and finally at 10~ 3 g cm -3 . As with the 
IBoss & Mvhill Il992t) initial conditions, the temperature differs by more 
than an order of magnitude from that given by the barotropic equation of 
state (dashed line) at late times. 



Figure 7. The evolution of the maximum temperature and den- 
sity during the calculation s that start from the initial conditions of 
IBoss & BodenheimeJ <1979l) . including the m = 2 density perturbation and 
rotation, for two different resolutions. The calculations were performed us- 
ing 50,000 (dotted red line) and 150,000 (solid green line) particles. The 
two resolutions give similar results. For comparison, the dashed line gives 
the temperature-density relation of the barotropic equation of state used by 
lBatdll998t) . 



ature versus density is much cooler than that from its spherically- 
symmetric counterpart (see Figure|5| red line). In fact, coinciden- 
tally, it is actually qu ite close to the evolution obtained with the 
IBoss & Mvhilll fl992) initial conditions. There are likely to be two 
reasons for the cooler temperatures at the same maximum density 
in the rotating collapse. First, the collapse to form the stars is de- 
layed by the rotational support (they form at t = 1.30 initial cloud 
free-fall times instead of just over one free-fall time for the spher- 
ically symmetric case). The slower collapse allows more time for 
energy to be radiated away from the object. Second, the presence of 
rotation means that the distribution of material around each of the 
collapsed objects is no longer spherically symmetric. In particular, 
there are low density cavities along the rotation axes and high den- 
sities in the discs. Radiation can more easily escape through these 
cavities than in the spherically symmetric case. 

We also teste d the effect of resolution on the 
IBoss & Bodenhe imeJ 1 19791) binary star initial conditions, 
using 50,000 and 150,000 particles (Figu re |7] red and green lines 
respectively). As with the Boss & Mvhilll iT992) initial conditions, 
increasing the resolution above 50,000 particles has little effect 
on the temperature evolution showing that the calculations are 
essentially converged. 

The difference in evolution between the radiation and the gas 
temperatures is shown in Figure |S| in a manner similar to Figure 
|2| Again the only significant difference occurs during the transition 
from the optically thin to optically thick regimes at densities around 
gem" 
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Figure 8. The evolution of the maximum gas (solid green line) and ra- 
diation (dotted red lin e) tem peratures at the maximum density during the 
IBoss & BodenheimeJ 1 19791) calculation, including the initial m = 1 den- 
sity perturbation and rotation. The calculation was performed using 150,000 
particles. The gas and radiation temperatures are identical except during 
the transition from optically thin to optically thick. For comparison, the 
dashed line gives the temperature-density relation of the barotropic equa- 
tion of state used by Bat^ |l9 



4 CONCLUSIONS 

The results of our study show that the initial conditions in a molec- 
ular cloud core such as the density and velocity configurations have 
a large effect on the evolution of the temperature as the cloud col- 
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Figure 9. Snapshots of the logarithmic column d ensity, N, (upper panels) and mass-weighted gas temperature, T g , (lower panels) parallel to the rotation 
axis for the 150,000 particle lBoss & Bodenheimej] ] 1979ft calculation, including the initial m = 2 density perturbation and rotation. The density scale covers 
-1 < log N < 3.5 with N measured in g cirT 2 while the temperature scale ranges from 1 < log T g < 3. The snapshots are given when the maximum densities 
are 10~ 12 , lO" 10 , and 10" 1 g cm . The corresponding times are 1.18, 1.24 and 1.30 initial cloud free-fall times. The molecular cloud core collapses to form a 
protostellar binary system. Note that the dense gas within the spiral arms is cooler than the surfaces of arms (e.g., lower middle panel and the left-hand object 
in the lower right panel). In the lower right panel, the protostar on the right has formed a stellar core and thus heated the surrounding gas significantly more 
than the protostar on the left. At the time plottted, the left protostar has not quite undergone molecular hydrogen dissociation to form a stellar core. The two 
objects form stellar cores at slightly different times because of small asymmetries in the initial conditions. 



lapses. The evolution of the maximum temperature and density 
during the collapse cannot be approximated by a single barotropic 
equation of state. Furthermore, even if a satisfactory description of 
this evolution was able to be formulated by a simple equation of 
state, it would do a very poor job of setting the temperature away 
from the location of the density maximum. We find the barot ropic 
equat ion of state used in many previous simulations (e.g. iBata 
1 19981) may underestimate the temperature by up to a factor of two 
to three for the density maximum, and by more than an order of 
magnitude in other parts of the cloud. In principle, this may seri- 
ously affect the evolution of a molecular cloud core, even altering 
how it fragments since the Jeans mass scales with gas temperature 
as Tg /2 . Therefore, if possible, radiative transfer should be used. 
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